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ABSTRACT 



We use CFHT/MegaCam data to search for outer halo star clusters in M33 
as part of the Pan- Andromeda Archaeological Survey (PAndAS). This work ex- 
tends previous studies out to a projected radius of 50 kpc and covers over 40 
square degrees. We find only one new unambiguous star cluster in addition to 
the five previously known in the M33 outer halo (10 kpc < r < 50 kpc). Al- 
though we identify 2440 cluster candidates of various degrees of confidence from 
our objective image search procedure, almost all of these are likely background 
contaminants, mostly faint unresolved galaxies. We measure the luminosity, color 
and structural parameters of the new cluster in addition to the five previously- 
known outer halo clusters. At a projected radius of 22 kpc, the new cluster is 
slightly smaller, fainter and redder than all but one of the other outer halo clus- 
ters, and has g' ~ 19.9, {g' — i') ~ 0.6, concentration parameter c ~ 1.0, a core 
radius Vc ~ 3.5 pc, and a half-light radius ~ 5.5 pc. For M33 to have so 
few outer halo clusters compared to M31 suggests either tidal stripping of M33's 
outer halo clusters by M31, or a very different, much calmer accretion history of 
M33. 



Subject headings: galaxies:individual (M33) 



1. Introduction and Background 



Globular cluster systems (GCSs) are important tracers of galaxy formation and 
evolution. For example, the substructure within a galactic halo reveals its merger 
history, and globular clusters (GC s) can be used as one tracer of such substructure (e.g.. 



Lynden-Bell fc Lynden-Bell 



19951 ). In this paper, we look at the Triangulum Galaxy (M33) 
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which is the third most massive galaxy within our Local Group although it is much less 
well-studied than the Milky Way (MW), the Magellanic Clouds or Andromeda (M31). 

In addition to the well-known stream from the disrupting Sagittarius dwarf galaxy, 
other evidence for substructure within the Milky Way in cludes the Monoc e ros ring, the 



Orphan stream, and other more su btle overdensities (e.g. 



Ibata et al. 



1995 



2002 



Belokurov et al. 



2006 



Newberg et al. 



20071 ). M31's substructure is being revealed in more and 



more detail with large-sc ale structures o : 



.1 

arcs, shells and strearns f|Ferguson et al. 



Richardson et al 



2008 



very 



McConnachie et al. 



low surface brightness 



20021: lib? 
aL 20oX 



Ibata et al. 



2005 



2007 



including several 



Kalirai et al. 



2006 



Subdivisions ^yithin the MW GCS have been observed (e.g., ISearle Sz Zinrull978 



Mackey fc Gilmore 



2004j ) with evidence that at l east some are the result of accretions 



of dw arf satellite galaxies (IBellazzini et al. 



2003 



Mackey &: Gilmore 



2004 



Forbes et al 



20041 ). Certain clusters still appear to be associated with their accreted satellit es 



most prominen 


tly, clusters associatec 


. wit : 


Newbere et al. 


2003; 


Bellazzini et al. 


2003 



AM-1, first discovered by 



1 Sagittarius (e.g.. 



Layden fc Sarajedini 



2000 : 



2003 . The most distant GC known in the MW is still 



Madore fc ArpI (119791 ) at a galactocentric distance of ^120 kpc. 



In other galaxies, clusters with large galactocentric radii (~ 120 kpc) reside in the My - r/j 



parameter space between 
now well-established (e.g., 



alomar-type clusters and ultra- 



Huxor et al. 



2005 



Gomez et al. 



aint dwa rfs, and this over 



2006 



and 



Belokurov et al. 



ap is 



20071). 



Within M31, there are now over 60 known clusters with a proiected radius greater 



than 30 kpc ( 



Mackey et al. 



iuxor et al. 



2005 



Martin et al 



2006 



Mackey et al. 



2007 



Huxor et al 



2008 



20101 ). Some of these distant clusters are rather unl ike their MW 



count erparts 



)3). 



20071). Most 



as they are found to be both more luminous and have larger sizes (iMackey et al.l 
recently it has been shown that the outer halo clusters appear to follow other substructure 
(streams of enhanced surface brightness), with the probability of chance alignment less than 
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1% (IMackey et al.l 120101 ) . iMackey et al.l conclude that the maiorit y of these clusters ar e 



accreted along with their host satellite galaxy, as first proposed by 



Searle fc Zinnl fll978h . 



Observations of the M33 clusters - bo t h you ng and old - have been collated 



in the catalogue by 



Sarajedini fc Mancond ( 120071 : SM hereafter). This catalogu e 



includes cluster identifications and data from ground-based observations ( 



Melnick fc D'Odorico 



1978 



Christian &: Schommer 



1982 



1988 



Hiltnei 



Mocheiska et al 



1960 



1998h. HST 



imag i ng ( 



2007 



Chandar et al 



Stonkute et al. 



1999 



2008 



identified clusters from 



2001 



Bedin et al 



Huxor et al 



Ma et al. 



(12001 



20091: 



2005 



Park fc 



San Roman et al. 



2002a 



I 



2007 



Sarajedini et al 



20091 ). and further data on 



2004al Jbl). The SM catalogue contains 



595 objects of which 428 are classified as high- confidence clusters (based on HST and 
high-resolution ground-based imaging). The most recent work, curre ntly not within th e 



SM ca talog ue, includes work based on CFHT/MegaCam i maging by 



(|2008|) and 



San Roman et al. 



(120101 ) and HST imaging by 



Zloczewski 



Zloczewski &: Kaluznyl (120091 ) 



et al. 



that contain 3554, 599 and 91 new star cluster c andidates, respectively. 



(All of these M33 



studies cover only the inner one square degree. 
~20% of the 3554 cluster candidates identified in 



Zloczewski fc Kaluznvi (20091 ) claim that 



Zloczewski et al. 



(120081 ) are hkely to be 



genuine clusters. Unlike 



the GCSs of the MW and M31, M33 is host to intermediate-age 



clusters ( ISarajedini et al. 



1998 



Chandar et al. 



20021 ). suggesting that the evolution of M33 



was different from that of both the MW or M31. 

Studying the Local Group gives us the best chance to observe the remnants of galaxy 
formation in detail, but M33 remains to be scrutinized in as much detail as either of its 
larger neighboring galaxies, or the Magellanic Clouds. The work on the M33 GCS has 
so far been constraine d to the class i cal dis k regions, with the exception of the four outer 



halo clusters found bv iHuxor et a. 



and one cluster by 



Stonkute et al 



(|2009[) between projected radii of 9.6 and 28.5 kpc 



(120081 ) at a projected radius of 12.5 kpc. The outer 
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halo clusters are import ant, not least because the most d istant clusters may be the last 



that were accreted (e.g. 



Mackey fc van den Bergh 



that M31's outer halo is rich with clusters. 



2005L 



Mackey et al. 



Huxor et al 



(120101 ) have shown 



(|2009[ ) undertook a search for 



M33 outer halo clusters through 12 sq. degrees of the Isaac Newton Telescope Wide-Field 
Camera data reaching to V~24.5 and i~23.5. The PAndAS data allow this search to be 
extended to larger radii, deeper depths and better image quality. This is the project that 
we undertake in this paper. We define outer h alo clusters to be thos e which are projected 



beyond the isophotal radius of M33 (~9 kpc, 



Cockcroft et al 



20111 ). Such objects are 



sufficiently remote that they are unlike ly to be associated with the main disk component of 



the galaxy: 



McConnachie et al. 



(120101 ) find little evidence from direct stellar photometry 



that the disk extends beyond that point. For co mparison, the isopho tal radius of NGC 253 



20031 ). Ultimately however, 



an Sc-type galaxy of similar size, is r ~ 9.8 kpc (iJarrett et al.l 
we will require metallicity and velocity measurements to determine more definitely whether 
these clusters belong dynamically to the disk or halo. 



2. Observations, Data Reduction and Calibration 



We use 41 im ages each in q' and i! tha t are part of the Pan- Andromeda Archaeological 



Survey (PAndAS; 



McConnachie et al 



20091 ) and were taken with the Canada-France-Hawaii 



Telescope (CFHT) /MegaCam which has a one square degree field-of-view. PAndAS includes 
over 300 sq. degrees, covering a region of sky that extends to a projected radius of 150 and 
50 kpc around the Andromeda (M31) and Triangulum (M33) Galaxies, respectively. 

Each of the 82 processed fields, themselves a stack of four or five raw images, is labelled 
according to the notation shown in Figure [H Fields M72 to M76, the five located along a 



line towards M31, were taken first (llbata et al 



20071 ): we note that for these, CCD chip 4 



was only working for M72, but not M73 to M76. The other fields M3301 to M3335 were 
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taken in subsequent runs. T 
archived images prepared by 



le central i r aage 



Ibata et al. 



(called M33c) was a composite of CADC 



fl2007h 



All M33 images were taken with sub-arcsecond seeing in both g' and i' . The average 
seeing on g' frames was 0.75" (standard deviation of 0.11"), and 0.66" on i' (standard 
deviation of 0.13"). Images have a resolution of 0.187" /pixel, and limiting magnitudes 
of g' ~ 25.5, i' ~ 24.5 ( AB mags on the SPSS sca le) at an S/N = 10. These data were 



previously presented in 



McConnachie et al. 



(|2010| ) who studied the stellar structure of 



the outer regions of M33 and found a large substructure, shown in Figure [H The data 
were pre-processed with Elixii^ by the CFHT team, and then reduced at th e Cambridge 



Astro nomical Survey Unit through a pipeline adapted for MegaCam images ( llrwin fc Lewis 
200l|). 



3. Cluster Search Methods 



We used two distinct methods to search for clusters within the M33 images: an 
automatic search and a separate follow-up visual inspection. In both cases, we started with 
searching the images that form an annulus in the middle of the frames around M33 (i.e., 
frames M3306 to M3316, including M74), followed by the images outside this annulus and 
finally the innermost frames. We chose this order so that we started in reg ions where the 



Huxor et al 



20091), leaving the 



crowding was low but evidence for clusters existed (e.g. 
inn ermost crowded fields until la st. We identified and marked all the confirmed clusters in 
the ISarajedini fc Mancond ( 120071 ) catalogue on the PAndAS images, to gain experience of 



their appearance. 



^ http : / / www . cfht . hawaii . edu/ Instrument s / Imaging / MegaPrime /dat apro cessing . html 
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3.1. Automatic Search 



We used Source Extractor (SE; iBertin &: Arnoutall996l ) to identify all objects in both 



g' and i' frames. Some values within the configuration files were changed to optimize 
finding clusters whilst cutting out as much contamination as possible (see Tabled]). We 
then converted from a chip-oriented pixel coordinate system to a world coordinate system 
for each object, using the IRAF/wcsctran routine, before matching the objects across the 
g' and i' frames using the IRAF/xyxymatch routine. As the latter routine only outputs 
object coordinates, we then re-assigned all SE parameters to the matched objects, so that 
we could apply selection criteria using a combination of magnitude, color, half-light radius, 
and ellipticity as measured by SE to pick out the cluster candidates. After numerous initial 
tests and iterations, we have adopted the following set of criteria: 

10.5 <g' < 14.5, (1) 



1.1 <^'-i' < -0.175* ^' + 3.4375, (2) 



e < 0.375, (3) 



3.5 < Tfi^, < 16.0, (4) 



and 



Tfiu. < -2.125* ^' + 41.5, (5) 
where e is the ellipticity(l - minor/major), g' and i' are the automatic magnitude values 



- 10 - 



returned from SEo, and tjiux is the radius (in pixels) estimated to enclose half the flux. The 
cluster candidates we select satisfy all five of the criteria. 



These selection criteria are shown in Figure El Th e 



they included almost all of the 



boun dary lines were chosen so that 



Sarajedini fc Mancond (120071 ) catalogue confirmed clusters 



at the edge of the disk (in the MegaCam fields M3301 and M76), while cutting out most 
of the contaminating objects such as stars and background galaxies. The central field, 
however, provides a special challenge because of the complex structure of the background 
light and differential reddeni ng. As can be seen in Figure [Sj our parameter boundaries do 



Sarajedini &: Mancond high confidence clusters (see Section [T]) 



not include every one of the 
in the central field. 

However, our aims here were specifically to isolate candidate clusters in the halo 
regions. Other types of objects (especially background galaxies) populate all areas of the 
three parametric diagrams in Figures [2] and [3l and after many iterations we adopted the 
boundary lines shown as a compromise between excluding contaminants and including real 
clusters. Nevertheless, the unavoidable fact is that our survey area is so large (more than 
40 square degrees around M33) that even our most careful objective search criteria leave a 
very large number of field contaminants, which dominate the numbers of objects found in 
the range of magnitudes, colors, and sizes that we are looking for. 

We produced a small thumbnail display region to identify each object selected by the 
above criteria. These regions were then displayed in the g' frame. Each object within these 
regions was then inspected visually and classified following the description in Section 13.21 
with 1 (high confidence cluster), 2 (possible cluster), 3 (background galaxy), 4 (unknown 



^The following are approximate conversions between true color-corrected magnitudes and 
SE automatic magitudes: g^^^^ = g'g^ + 6.2 and zj^^g = i'g^ + 6.4. 
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object) or 5 (stellar object). Examples of objects in the categories 1, 2 and 3 are shown in 
Figure HI 



3.2. Visual Inspection 



The next stage in our classification procedure was, following Section 13. H to inspect all 
objects that had not already been selected by the automated criteria, i.e., we looked at 
all SE detections that did not fall within the se l ection boxe s shown in Figure [2J This was 
similar to the method employed by 



Huxor et al. 



(120081 



20091 ). Only i' frames were inspected 



this way; red-giant branch stars are brighter in i' than in g' and so clusters, and RGB stars 
in their halo, would appear more obvious by their resolution into stars (no background 
galaxies would be resolved into individual stars). By also conducting a visual search in 
addition to the search via the selection criteria, we ensured that any obvious cluster or 
candidate cluster would not be overlooked, and also ensured that all g' and i' frames will 
have been inspected. 

The easiest objects to classify were the obvious clusters and background galaxies. 
Clearly-resolved clusters appeared as having a circular or slightly elliptical core, with 
uneven contours and resolved stars around the central core. Less obvious were group 1 
objects which had slightly uneven contours. The least obvious candidates, group 2 objects, 
were the compact objects that could be clusters or galaxies, and as a result the numbers 
of "possible clusters" were the greatest especially in the central regions. As the contrast 
and scaling were changed, some objects smoothly grew, some were a faint smear with no 
sharp edges, and others could be seen to display spiral shape. If the object had smooth 
contours and if it was in a group of other objects that were clearly galaxies, the object was 
likely to be a galaxy and not a star cluster. Group 4 objects did not look like a cluster, a 
cluster-candidate, a galaxy or a star. 
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As noted above, the influence of background contamination by galaxies on tliis 
selection process should not be underestimated. In essence, this is a needles-in-a-haystack 
process where we are attempting to find a small number of clusters in a huge population of 
contaminants, and even though our selection and culling is rigorous, there remain a large 
number of objects whose nature is ambiguous from the current data. Higher resolution 
imaging, imaging in the near infrared where the cluster red giants would be better resolved 
(and which also can have better seeing), or ultimately spectroscopy, will be required for 
more definitive elimination of the last contaminants. 



Results 



There was only one definite new outer halo cluster discovered in our study at a 
projected radius of 87" (or 22 kpc, assuming a distance to M33 of 870 kpc). It was 
found using the auto mated search . The new cluster is named M33E following the naming 

i ^our o f the five previously-known outer halo 



conventi on begun in 



Huxor et al. 



clusters ( Stonkute et al. 



2008 



(120091). 



Huxor et al. 



20091 ) were easily recovered. Cluster D was 



identified but was too compact to have been recovered without pri or knowledge. C 



to E and S are shown in Figure where S is the cluster found by 



Stonkute et al. 



usters A 



(120081 ). 



There were 2440 candidates spread throughout the M33 halo; that is, in the region 
outside of the central MegaCam image. 87 (5) highest-confidence cluster candidates and 
2294 (54) possible clusters were found by the automated (visual) search method. 

The numbers of all classified objects from both the automated and visual inspection 
searches are shown in Table El Results of the above searches were plotted within the original 
selection criteria, and are shown in Figure [61 We wanted to exclude the maximum amount 
of parameter space so that we could increase the efficiency of the automated search, and it is 
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not obvious from this figure that more space could have been excluded. Radial density plots 
for the categories 1 (high confidence cluster), 2 (possible cluster), 3 (background galaxy) 
are also shown in Figure [71 We compared these number densities at large radii to control 
fields, also taken with MegaCam and with very similar image quality, from the M31 outer 
halo and the field near the Draco dwarf spheroidal. The M31 fields are two square-degree 
fields selected directly from the PAndAS data, at a similar Galactic latitude to M33 of 
-31.33 degrees, at the edge of the PAndAS footprint around M31 (i.e., at a projected 
radius of ~150 kpc) and did not contain any clusters - either previously-known clusters, 
or clusters detected in the PAndAS images. The Draco fields are seven square-degree 



fields at a Galactic latitude of 34.72 degrees (ISegall et al. 



20071 ). Our searches were again 



applied to the control fields following exactly the same selection criteria, and we obtained 
an average density of each category of objects in the control fields. The radial distribution 
plots indicate that few if any of the category 3 objects are genuine clusters since they show 
little detectable central concentration to the galaxy outside the crowded disk region. For 
all three categories plotted, the number density settles down to a virtually constant level 
similar to that of the M31 control fields for r > 1 degree, consistent with the conclusion 
that there are few clusters left to be found in the M33 halo down to the PAndAS limiting 
magnitudes. (The number density of all objects in the Draco fields is much lower than that 
in either the M33 or M31 fields, highlighting that it was appropriate to compare M33 with 
the M31 control fields.) If we count the number of candidates for the combined objects of 
classes 1 and 2 in Figure [7] for r > 10 kpc, and then subtract off the M31 background, we 
are left with approximately 210±130 candidates (the error is estimated using the error on 
the M31 background). This number is simply an estimate of the outer halo clusters that 
possibly remain to be discovered, using the data we have in hand. 210 clusters would be a 
generous upper limit, given the field contamination issues that we discuss. 

We next measured the g' and i' magnitudes, and the colors of the six outer halo 
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clusters. The results are shown in Table |3l For clusters A, B, C and S we use an aperture 
radius of 40 pixels (7.5") and a sky annulus between 60 and 80 pixels (11.2"-15.0"); 
for the smaller clusters we used 20 pixels (3.7") with a sky annulus between 20 and 40 
pixels (3.7"-7.5") for D, and 30 pixels (5.6") with a sky annulus between 50 and 70 pixels 
(9.4"-13.1") f or E. We assurne an extinction correction of 0.16 in g' ^ 0.09 in z', 0.14 in 1/, 



and 0.08 in / fISch 
and 



Huxor et al 



egel et al. 



19981 ) ■ and a distance of 870 kpc to M33, consistent with SM07 



(120091 ). and corresponding to a distance modulus of (m-M)o = 24.69. We 



note that th ere is some disagreement in the literature regarding the distance to M33; see 



references in 



McConnachie et al 



( j2010f ). For the magnitude and color conversions from 



((?', i') to iy, I) we used 



V = g- (0.587 ± 0.022)(^ - r) - (0.011 ± 0.013) 



(6) 



and 



I = i- (0.337 ± 0.191)(r - i) - (0.370 ± 0.041) 



(7) 



from 



Chonis fc Gaskelll (120081 ). and 



(r-z)o = 0.37(c/ - r)o + 0.006 



from iBihr et all toO& l . 



Compar i ng th e magnitudes and colors that we measure in this paper with those in 



Huxor et al. 



(120091 ) for clusters A, B, C and D, we find differences of Vo< 0.3mag and 



(y — /)o < O.lmag. For cluster S, 



Stonkute et al. 



(120081 ) measures a V magnitude of 18.5 



(at ~ 2r/i, which roughly corresponds to our annulus size). As we measure Vo~18.5, the 
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difference between our measurements is the extinction value of 0.14. 



Our crude estimate for tlie cluster magnitude limit is currently gl-^ ^ 20 {Mg ^ —4.8). 
We will quantify this more accurately in an upcoming paper by inserting fake clusters 
and testing recovery rates using our search methods. Although our current search limit is 
comfortably faint, there are small numbers of still less luminous clusters known to exist in 
the Milky Way, for example (the faintest, sparse Palomar-type objects; see Figure |9]). We 
can therefore place no quanti tative limits on th e numbers of such objects yet to be found 



in M33. Note that cluster D ( IHuxor et al. 



20091 ) is a magnitude fainter than our estimated 



limit, but was found with HST imaging. We would not expect to recover such a cluster 
independently with the MegaCam data. 

Finally, we measured the structural parameters of all six outer clusters, including 



the concentration parameter, and co r e, ha. 



code described by 



McLaughlin et al 



f-light and tidal radii. We use the GRIDFIT 



( 120081 ) ■ which fits various King-type cluster mode ls 



convolved with the measured 



Kind (119661) and 



Larsen et al. 



(l2002h with the 



'SF to each object. Here we attempt to fit 



KingI (119621) 



WilsonI (Il975[) models t o each object. We also use the KFIT2D code of 



King! (119661 ) model as an independent measure. The results 
of all fits are shown in Table HJ and examples of the fits are shown in Figures [HI We also 
include an independent measurement of the half-light radii using the curve of growth of 
the clusters (rap). For A and D, not all of the three models converged to successful fits, 
but the other four clusters gave high consistency among themselves for their radii. Fitting 
models to cluster A was not successful because of its diffuse nature, while cluster D was 



extremely sma. 



Stonkute et al. 



1. We also note that cluster S has an unusual feature in its surface profile; 



( 120081 ) also note that this cluster is asymmetrical in its in ner regions. 
Comp aring our measured structural parameters for cluster S with those from iStonkute et al. 
(120081 ) we similarly find that this cluster has a very large core radius, although we measure 
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slightly smaller quantities for each radius - even adjusting the value for our assumptions for 
the distance to M33. 

In Figure m we show the locations of all six M33 halo clusters in luminosity versus Vh, 
compared with the Milky Way GCs. All clusters have low concentrations, similar to the 
Palomar outer halo clusters in the Milky Way. Their half-light radii range from 4 to 20 pc, 
all larger than the typical mean r/^ ~ 3 pc for the standard Milky Way clusters, but placing 
them in a similar range as many of the outer halo Milky Way clusters. A full comparison of 
all M33 clusters will be done in a future paper. 



Discussion and Conclusions 



We search more than 40 square degrees of the halo of M33 with CFHT/MegaCam data 
for outer halo clusters using both an automated search and visual inspection. Unexpectedly, 
we find only one new cluster, which is sm aller, fainter and slightly redder than the three 



INT clusters found by iHuxor et al. 



(120091 ) a nd the cluster found by iStonkute et all (120081 ). 



However, it does lie within the INT area of 



Huxor et al. 



(120091 ) but the object was not 



previously recognized due to its small size and faint luminosity. At a projected radius of 22 
kpc, the new cluster has g' ~ 19.9, {g' — i') ~ 0.6, concentration parameter c ~ 1.0, a core 
radius ~ 3.5 pc, and a half-light radius ^ 5.5 pc. Its proj ected location is close to the 
feature observed in the stellar substructure (see Figure [T] and 



McConnachie et al. 



Huxor et al. 



20101). 



(120091 ) note that the mean color of the previously-known outer halo clusters is 
slightly redder {{V — I)o = 0.88 ± 0.05 mag) than the inner clusters {{V — /)o = 0.67 ± 
0.30 mag). Our new cluster is redder still by ~0.2 mag. 

M33 has only six definite outer halo clusters between projected radii of 9 kpc < r < 50 
kpc and to g'^^ ~ 20. We also find 2440 cluster candidates of various degrees of confidence. 
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and although the vast majority are hkely to be background contaminants, at least some 
of the ~ 90 highest-confidence candidate objects beyond the M33 disk may be faint but 
genuine clusters. We cannot yet assume all the highest-confidence candidate objects are 
clusters without further confirmation. We will use IR data (now being acquired) and 
structural parameters in an upcoming paper to determine this more securely. 

How many clusters could we expect to find in M33? M31 has 67 outer halo clusters 
already discovered, 61 of which lie in the PAndAS footprint that has been analyzed so 
far. These clusters have comparable luminosity to the M33 outer halo clusters and are 



located at projecte d radius 30 kpc < r < 130 kpc (IMackey et al 



2010 



Huxor et al 



Huxor et al. 



201 ih . 



(120091 ) found a GC surface density of ~0.4 deg~^ with their 12 deg^ study, 
which they note is about half that derived for M31 over the radial range 30kpc < r < 130 
kpc. Here we find an even lower GC surface density of 0.15 deg~^. We note that the search 
in M3rs outer halo is not yet complete so its GC surface density is likely to increase. M33 
appears to therefore lack this type of cluster. We briefiy mention two scenarios that could 
have resulted in this observed difference. 

M33 could have had a different accretion history co mpared to M31 - a conclu sion that 



San Roman et al. 



2010|), but is 



has been drawn before from studies of the inner regions (jS 
now also indicated by the outer halo data. If M33 never interacted with M31 before, M33 
would have had a dramatically less active accretion history. 



The most compelling evidence for a n accretion 



from the Sagittarius d warf in the MW ([Ibata et al 



origin for the outer halo clusters comes 



19951 ). and from the GCs and tidal 



20101 ) , but it is still far from clear how general a result 



debris streams in M31 (IMackey et al. 
this is." 

However, another exciting and more likely prospect, given the tidal distortion of M33, 
is that perhaps some of M33's outer halo clusters were heavily stripped off in a previous 
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dynamical interaction with M31 (iHuxor et al 



2009 



San Roman et al. 



2010h . Some of the 



GCs originally belonging to M33 may now be closer to M31, but it will be difficult to 
disentangle the populations. A more detailed comparison will require spectroscopic studies 
of these clusters to determine properties that may link the divided populations. Although 
unlikely, some clusters may be beyond the area that we have imaged so far around M33. 
These scenarios are not mutually exclusive. Further discussion and comparison with the 
M33 halo star population will come in subsequent work now in progress. 
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Table 1. Source Extractor values used on the MegaCam images. 



Parameter 


Default Value 


New Value 


DETECT3IINAREA 


5 


4 


THRESH.TYPE 




RELATIVE 


DETECT.THRESH 


1.5 


5 


DEBLENDJNTTHRESH 


32 


8 


CLEAN.PARAM 


1.0 


1.5 


PHOT.APERTURES 


5 


3 


PHOT.AUTOPARAMS 


2.5,3.5 


2.0,2.5 


PHOT_PETROPARAMS 


2.0,3.5 


2.0,2.5 


PHOT_FLUXFRAC 




0.5 


SATUR.LEVEL 


50000 


60000 


MAG.ZEROPOINT 


0.0 


(/'=26.7,i'=25.98 


BACKPHOTO.TYPE 


GLOBAL 


LOCAL 
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Table 2. C ategorized obiects in M 3 3 PA ndAS frames. The column headers are as follows: 



SMI are the lSarajedini fc Mancond (j2007l ) catalogue's confirmed clusters, OH are the outer 
halo clusters, 1 are our highest- confidence clusters, 2 are the possible clusters, 3 are the 
background galaxies, 4 are unknown objects and 5 are stellar objects. The numbers in 

brackets after the highest-confidence and possible clusters indicate those candidates which 

matched objects in SMI. 



Frame SMI OH 1 2 3 4 5 



M3301 


34 


- 


7(3) 


354(16) 


106 


114 





M3302 


- 


- 


1 


69 


56 


8 


27 


M3303 


1 




3 


96 


43 


16 


8 


M3304 


- 




2 


69 


86 


16 


7 


M3305 


- 


j^6, e 


2 


83 


51 


19 


5 


M3306 


- 


- 


2 


49 


92 


20 





M3307 


- 




2 


31 


79 


16 





M3308 


- 


- 





45 


149 


13 





M3309 


- 


- 





77 


104 


11 





M3310 


- 


- 





64 


77 


7 





M3311 


- 


- 


1 


27 


74 


6 





M3312 






3 


38 


96 


11 





M3313 






2 


64 


60 


13 





M3314 






3 


44 


55 


13 


1 


M3315 






2 


35 


87 


31 





M3316 









56 


81 


14 





M3317 






1 


44 


114 


9 





M3318 









25 


121 


5 





M3319 









55 


96 


16 





M3320 









29 


86 


6 





M3321 









8 


74 


9 





M3322 






3 


26 


66 


11 


1 


M3323 









18 


69 


6 


1 


M3324 






1 


63 


61 


16 


3 


M3325 









31 


61 


15 


3 


M3326 






1 


29 


104 


15 





M3327 









24 


81 


14 


3 


M3328 






1 


32 


115 


14 


3 


M3329 






5 


72 


81 


19 


6 


M3330 






5 


67 


117 


15 





M3331 






2 


61 


127 


25 





M3332 






1 


58 


148 


22 


6 
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Table 2 — Continued 



Frame 


SMI 


OH 1 


2 


3 


4 


M3333 




3 


39 


129 


16 


M3334 




2 


35 


116 


15 


M3335 




2 


29 


87 


25 


M72 




2 


28 


71 


11 


M73 




1 


45 


93 


13 


M74 







52 


61 


12 


M75 







43 


95 


22 


M76 


19 


2"''= 32(8) 


234(2) 


54 


86 


M33c 


374-f 


259(95) 


1521(138) 


84 


954 



Also found in M3305 
* Also found in M76 

Also found in M3304 

Also found in M3303 

Also found by Huxor et al. (_2009) 
f Does not include those on frames M3301, M3303 or M76) 
^ New cluster identified in this paper. 



Table 3. Outer halo cluster positions, lum inosities anc 



870 kpc to M33, consistent with SM07 and 



Huxor et al. 



colors. We assume a distance of 



fcoOOh . and that M33's center is 



located at (0l'*33™50.9^ 30'='39"^37 



Degrees 



Galactocentric distance 



Cluster 


RA 


Dec 


arcmins 


kpc 


9'o 


(g' - i> 




{V-I)o 


A 


23.92388 


28.82086 


112 


28.4 


19.1 


0.7 


18.8 


0.9 


B 


24.00865 


29.96372 


48 


12.2 


17.8 


0.7 


17.5 


0.8 


C 


24.31026 


31.07433 


45 


11.3 


18.4 


0.8 


18.1 


0.8 


D 


23.75916 


31.23925 


37 


9.4 


21.3 


0.8 


20.9 


1.0 


E 


23.84466 


32.07559 


87 


21.9 


19.8 


0.5 


19.6 


1.1 


S 


23.24374 


29.8675 


49 


12.3 


18.9 


0.8 


18.5 


0.8 
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Table 4. Outer halo cluster structural parameters, including concentration, c, core radii, 
half-light radii, and tidal radii, using both GRIDFIT and KFIT2D. We report the best-fit 
to the data, whether in g' or in i' . Also shown are half-light radii estimates (r^p) using 



curves-of-growth, for an independent check on the upper limit o 
assume a distance of 870 kpc to M33, consistent with SM07 and 



the half-light radi i. We 



Huxor et al. 



fl2009h . and 



corresponding to a distance modulus of (m-M)o = 24.69. 



Radii (pc) 

Cluster Model Band Seeing/FWHM c Core Half Tidal r(ap) 





K62 


i 


3.8 


0.5 


16.9 


20.1 


68.0 


11.7 


A 


K66 


i 


3.8 


0.4 


17.7 


20.3 


83.2 






W 


















kfit2d 


i 


3.8 


0.7 


11.7 


11.6 


59.8 






K62 


g 


4.0 


0.8 


6.1 


9.4 


44.0 


9.0 


B 


K66 


g 


4.0 


0.9 


6.3 


9.4 


56 






W 


g 


4.0 


1.0 


6.6 


9.4 


87.3 






kfit2d 


g 


4.0 


1.1 


6.0 


10.9 


78.5 






K62 


g 


4.0 


0.8 


5.6 


8.7 


40.7 


7.8 


C 


K66 


i 


3.3 


1.0 


4.5 


7.3 


48.1 


7.8 




W 


i 


3.3 


1.0 


4.9 


7.1 


70.8 






kfit2d 


i 


3.3 


1.3 


3.8 


8.8 


81.9 






K62 


g 


4.5 


0.3 


4.7 


4.8 


13.6 


4.7 


D 


K66 


















W 


















kfit2d 


g 


4.5 


0.4 


4.1 


3.7 


9.8 






K62 


i 


2.6 


0.8 


3.3 


5.2 


25.5 


5.1 


E 


K66 


i 


2.6 


0.9 


3.5 


5.2 


31.4 






W 


g 


3.1 


1.2 


3.6 


5.6 


65.8 


5.1 




kfit2d 


i 


2.6 


1.1 


2.7 


5.0 


32.4 






K62 


g 


4.6 


0.3 


15.8 


15.7 


42.7 


18.7 


S 


K66 


i 


3.4 


0.4 


14.8 


16.7 


67.31 


18.7 




W 


i 


3.4 


0.2 


15.8 


19.2 


121.4 






kfit2d 


i 


3.4 


0.7 


17.9 


18.7 


93.7 
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S (degrees) 



Fig. 1. — (Left) The 41 PAndAS frames by CFHT/MegaCam around M33 used in this pa- 
per. Each square represents the one square-degree field-of-view of MegaCam. M33c is the 
central field. The two circles represent projected radii of 10 kpc and 5 kpc centered on M33. 



Also s hown are the locations of the high- confidence clusters in the 



(120071 ) catalogue (the outer halo cluster found by iStonkute et al. 



Saraiedini &: Mancone 



(120081) be ing the only 



catalo gue's object outside the 10 kpc radius), the four outer halo clusters in 



Huxor et al. 



(120091 ) enclosed in one box and the newly discovered outer halo c luster enclosed in a double 



box. The current (October, 2010) online 



Sarajedini fc Mancond cat alogue is further subdi- 



vided to show the location of the 296 original clusters in red from 



Sarajedini &: Mancone 



( 2007 ) in addition to the subseq uently-discovered 32 c 



(120071 ) . 115 cluster s in cyan from 



ters in gr een from IS an Roman et al 



found by 



Zloczewski fc Kaluznyl ( l2009l ) and 



usters in black from 



Park fc Lee 



oczewski &: Kaluznyl ( 120091 ) . and 161 (115 new) clus- 



( 20091 ). There is much overlap between the clusters 



San Roman et al. 



halo clusters overlaid on the substructure map from Figure 13 in 



2009h. fRight) Th e oute r 



McConnachie et al. 



(I2010h . 
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Fig. 2. — Two halo fields and objects on which Section I3.1f s selection criteria were based. 
The black dots are the objects of all kinds detected by SE. The green dashed lines show 
the boundaries of the selection criteria (equations [1] to |5]), and the green squares enclose 



those points which were picked out 
the high-confidence clusters in the 



jy all five of the selecti o n crit eria. The red squares show 



Sarajedini fc Mancond (120071 ) catalogue. The following 



are approximate conversions between true color-corrected magnitudes and SE automatic 



magitudes: g[j.^^ = g'g^ + 6.2 and 



true 



''SE 



+ 6.4. 
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Centre 




10 15 20 

g mag auto 



Fig. 3. — Objects within the central field. As with Figure 121 the black dots show the SE 
detections. The green dashed lines show the selection criteria, and the green squares enclose 



those points which were picked out 
the high-confidence clusters in the 



3V all five of the selecti o n crit eria. The red squares show 



Sarajedini &: Mancond (120071 ) catalogue. The following 



are approximate conversions between true color-corrected magnitudes and SE automatic 
magitudes: g[^^^ = g'^^ + 6.2 and i'^^^^ 



I'sE + 6.4. 
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Fig. 4. — Six examples of classified objects in g' (top row) and i' (bottom). From left to 
right: two examples each of 1 (high-confidence clusters), 2 (possible cluster candidates), and 
3 (galaxies) . More detail is apparent when changing the scale and contrast in a DS9 window. 
Each box is 20" square, corresponding to about 84 pc square (at 870 kpc). 




Fig. 5.— M33 A, B, C, D, E and S (left to right) in g' (top) and i! (bottom). Each box is 
20" square, corresponding to about 84 pc square (at 870 kpc). 
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a 




-1 1 

^g— i) mag auto 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 

4 6 8 10 12 14 

flux radius 



10 11 12 13 14 15 

g mag auto 



Fig. 6. — The classified objects after visual inspection shown in relation to the automatic 
selection criteria of Section I3UI Class 1 represents high confidence candidates, class 2 are 
possible clusters, class 3 are background galaxies, class 4 are unknown objects, and class 5 
are stellar objects. To calibrate the "mag auto" values, we added the zeropoint term and 
corrected for airmass and color. 
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Fig. 7. — Radial densities of objects in circular annuli from categories 1 (high-confidence 
clusters), 2 (possible clusters) and 3 (background galaxies) along the top, and the sum of 
1 and 2, and also 1, 2 and 3 on the bottom. Error bars on the points are simple ri^/^ 
uncertainties. Also shown are the mean and 1-sigma errors of the same categories of objects 
in control fields from the M31 outer halo (solid line) and Draco (dashed line). The Draco 
number densities are so low, as mentioned in Section HI that in two cases the 1-sigma errors 
are larger than the mean and so cannot be plotted on a log scale. In these cases, only the 
mean and upper 1-sigma error are shown. 
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Fig. 8. — Examples of radial profiles for each of the six outer halo clusters. The solid points 
indicate the data, the fine is the bestfit, the open points show the profile for the PSF. Note 
that for clusters A and D, the fits did not converge to a simple King- type model solution 
adequately. 
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Fig. 9. — Absolute m agnitudes versus the half-light radius in parsecs for both the Milky Way 
clusters (jHarrislll996l ). shown in small boxes, and the six M33 outer halo clusters, shown as 
larger boxes. 



